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Abstract 

The problem of extrapolating the series in powers of small variables to the region of 
large variables is addressed. Such a problem is typical of quantum theory and statistical 
physics. A method of extrapolation is developed based on self-similar factor and root 
approximants, suggested earlier by the authors. It is shown that these approximants 
and their combinations can effectively extrapolate power series to the region of large 
variables, even up to infinity. Several examples from quantum and statistical mechanics 
are analysed, illustrating the approach. 



Keywords: Power series; summability; self-similar approximants; bound states; compu- 
tational methods in statistical physics. 



PACS numbers: 02.30.Lt, 02.30.Mv, 02.60.Gf, 03.65. Ge, 05.10-a 



1 



1 Introduction 



A common problem in physics is how to extrapolate the series, derived by means of pertur- 
bation theory in powers of a small variable, to large values of this variable [1]. The known 
methods of extrapolation are mainly numerical and rather complicated [1,2]. Recently, being 
based on the self-similar approximation theory [3-8] and employing the algebraic self-similar 
renormalization and self-similar bootstrap [9-11], we have derived new types of approximants 
allowing for the summation of power series, which are the self-similar exponential approxi- 
mants [11-14], self-similar root approximants [15-17], and self-similar factor approximants 
[18,19], which, for brevity, could be called superexponentials, superroots, and superfactors, 
respectively. 

For the purpose of extrapolation, superexponentials are useful for describing the processes 
with exponential characteristics, such as financial time series [20-23] or the development of 
ruptures [24]. Extrapolation with the help of superexponentials was employed for analyzing 
and predicting financial markets [25-29] as well as ruptures and fractures of materials [30-33] . 

The processes, characterized by power-law behaviour, require for their extrapolation 
either the superroots or superfactors. Actually, both these constructions are the realization 
of the same resummation procedure under slightly different initial assumptions. It is the 
aim of the present paper to consider the extrapolation of asymptotic series by means of 
superroots and superfactors, and also by constructing their combinations. We especially 
concentrate on the problem of extrapolation of series, valid for small values of variable, to 
the region, where this variable tends to infinity. We illustrate the consideration by several 
examples typical of quantum and statistical physics. 

2 Factor and Root Approximants 

Suppose we are looking for a function f(x) of a real variable x, which is a solution of a 
complicated physical problem. Let this problem be so much complicated that it can be 
treated only by means of perturbation theory resulting in an expansion 

f(x)~J2a n x n (x-0) (1) 

n 

in powers of the asymptotically small variable x. But assume that our final goal is to find the 
behaviour of function f(x) at very large x, such that x — > oo. This is exactly the extreme 
variant of the extrapolation problem: how to extrapolate the finite series 

f k {x) = j^a n x n , (2) 

n=0 

having sense solely for x — * 0, to the region, where x — > oo? 

An effective summation of the power series (2) can be done in the frame of the self- 
similar approximation theory [3-8], resulting in self-similar root [15-17] or factor [18,19] 
approximants. We shall not repeat here the derivation of these approximants, since all 
details of the derivation procedure can be found in our earlier papers, but we shall just 
present the resulting formulas that will be employed in the following sections. 
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If 2k terms of a perturbative series are available, its effective summation can be done by 
the factor approximant 

/2*(z) = ao flil + Axp, (3) 
i=i 

where 2k control parameters A, and can be obtained by the accuracy-through-order 
procedure, expanding Eq. (3) in powers of x up to order 2k and comparing the like-order 
coefficients with those of the perturbative expansion f2k(x). This is the way of defining the 
even-order factor approximants. The odd-order approximants can be constructed in several 
ways [19], the simplest of which is by leaving untouched the zero-order term, which yields 

k 

fZ k+1 {x) = a + <n x J](l + A x) ni ■ (4) 

i=i 

Here and in what follows, the labelling of the factor approximants is done according to the 
nontrivial number of perturbative terms required for the approximant construction. Thus, 
an f k (x) approximant requires the knowledge of terms up to order k. 

Instead of prescribing one of the terms in the small-x expansion (1), one can impose a 
restriction on the behaviour of the sought function at large x — > oo, provided this behaviour 
is known, where 

f(x) ~ x 13 (x — > oo) . (5) 
This restriction for the even approximant (3) yields the condition 

k 

i=l 

while for the odd approximant (4), one has 

k 

i + J2^ = P ■ 

i=i 

Generally, several such conditions could be invoked, if one would know several terms describ- 
ing the behaviour of f(x) at large x — > oo. 

A &;-order root approximant can be defined through k terms of the large- £ expansion 
[15-17]. Here we shall consider another way of defining root approximants, by expanding 
them in powers of x and comparing such expansions with the given perturbative series (1). 
Then an even-order root approximant 

fg\x) = a (. . . ((1 + A 1 xr + A 2 x 2 )" 2 + . . . + A k x k ) nk (6) 

requires the knowledge of 2k nontrivial terms of the small-x series (1). The upper index in 
approximant (6) labels different possible solutions for the control parameters Ai and rii, since 
the accuracy-through-order procedure for Eq. (6) is not unique. This is contrary to the way 
of defining the superroot control parameters from the large-x expansions, which is a uniquely 
defined procedure [34]. An odd-order root approximant, with the control parameters defined 
from the small-x expansion, reads as 

if+i = go + a ± x (. . . ((1 + A, x)^ + A 2 x 2 ) m + ... + A k x k ) nk . (7) 
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Note that the control parameters Aj, and rii in each of the approximants given by Eqs. 
(3), (4), (6), and (7) are of course different, and we use the same letters for their notation in 
order to avoid too cumbersome nomenclature. In the present paper, we shall also consider 
the hybrid approximants combining the forms of factor and root approximants. 

In those cases, when the subsequent self-similar approximants f£(x) display substantial 
oscillations, it proved effective to introduce the averaged form 

T k (x) = X>(s) f*(x) , (8) 
i=i 

where f*(x) are weighted with the probabilities 

Pi(x) = ' 1 '] , , (9) 
in which N is the number of available approximants and 

are the mapping multipliers [35], where we may set f£(x) = a + a±x. Averages (8) can be 
defined for both factor as well as root approximants. 

3 Hybrid Factor-Root Approximants 

Here we show the way of constructing factor approximants, root approximants, and their 
various combinations. We illustrate this on the case of an asymptotic expansion obtained 
from the simple function 

f(x) = - ln(l +x) . (11) 

Ob 

Expanding this function in powers of x gives series (1) with the coefficients 

(-l) n 

- = • (12 > 

Retaining in this expansion the terms up to 6-th order, allows us to define the factor ap- 
proximant 

f*(x) = (1 + A 1 x)^(l + A 2 x) n2 (l + A 3 x) n * , (13) 

whose control parameters are to be found from the accuracy-through-order procedure, which 
yields 

A 1 = 0.9767, A 2 = 0.6261, A 3 = 0.1830 , 
ni = -0.3503 , n 2 = -0.1935 , n 3 = -0.2009 . 
At large x, approximant (13) behaves as 

f(x) ~ A^A^Af x "i+ ra 2+«3 ^ _ ? ( 14 ) 
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where 

A^A^Af = 1.5528 , m + n 2 + n 3 = -0.7447 . 
The root approximant of the same order is 

f^\x) = {[(l + A lX r+A 2 x 2 ] n2 +A 3 x 3 } n3 . (15) 

Expanding approximant (15) in powers of x, comparing this expansion with series (1), and 
equating the coefficients at like powers, we get the system of equations possessing 28 so- 
lutions. However, a natural restriction, limiting the number of admissible solutions, is the 
requirement that approximant (15) be real. Then only 5 solutions remain. The upper index 
j in Eq. (15) enumerates these solutions, which we shall analyse in turn. 
For the solution (x) we have 

A x = 0.9059 , A 2 = 0.0728 , A 3 = 0.0932 , 

ni = 2.2152 , n 2 = 0.9424 , n 3 = -0.2644 . 

At large x, we find 

ff\x) ~ Af x 3m (x -> oo) , (16) 

with 

Af = 1.8728 , 3n 3 = -0.7931 . 
For another solution fi 2 \x), the parameters are 

A t = 2.0426 , A 2 = 1.2353 , A 3 = 0.1938 , 

ni = 1.0002 , n 2 = 1.0031 , n 3 = -0.2440 . 
The asymptotic behaviour at large x is the same as in Eq. (16) but with 

Af = 1.4924 , 3n 3 = -0.7320 . 

The next solution fi 3 \x) has the parameters 

A x = 1.2509 , A 2 = 0.2772 , A 3 = 0.0135 , 

ni = 1.0616 , n 2 = 1.7956 , n 3 = -0.2097 . 
The asymptotic behaviour at large x is 

f!?\x) ~ A™ 2 " 3 x 2n2ns (x -> oo) , (17) 

where 

A n 2 « 3 = L6211 ; 2n 2 n 3 = -0.7531 . 
The parameters for (x) are 

Ai = 1.1673 , A 2 = 0.1682 , A 3 = -0.0545 , 
ni = 0.8629 , n 2 = 1.5029 , n 3 = -0.3303 . 
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The asymptotic form at large x is the same as in Eq. (17), but with 

A mn 3 = 2A227 ^ 2nsn3 = _ .9928 . 

For the approximant (x) , we have 

A x = 1.0761 , A 2 = 0.2497 , A 3 = 0.2389 , 

ni = 1.9117 , n 2 = 1.0210 , n 3 = -0.2381 . 
The behaviour at large x is of type (16), but with 

A3 3 = 1.4062 , 3n 3 = -0.7143 . 

All approximants f^\x) are close to each other. 

The hybrid approximants of the same order can be defined as follows. A possible form is 

$\x) = [(1 + Ai x) n ^ + A 2 x 2 }" 2 (1 + A 3 x) n * . (18) 

Defining the parameters Ai and rij by the accuracy-through-order procedure, we again meet 
the problem of multiple solutions. And again we reject those solutions that result in complex- 
valued or divergent functions (18), since the sought solutions must be finite and real for all 
< x < 00. The remaining approximants are again close to each other. Thus, for fjf\x) 
we have 

Ax = 0.8053 , A 2 = 0.1131 , A 3 = 0.9767 , 
n x = 0.9966 , n 2 = -0.1970 , n 3 = -0.3501 . 
The large-x behaviour is 

ff\x) ~ A 2 2 A 3 3 x 2n2+ns (x -> 00) , (19) 

where 

AVfAf = 1.5490 , 2n 2 + n 3 = -0.7441 . 
One more solution gives fQ 7 \x) of type (18), but with the parameters 

Ai = 1.2532 , A 2 = 0.2303 , A 3 = 0.7320 , 

ni = 1.0251 , n 2 = -0.2425 , n 3 = -0.2574 . 
At large x, one has the same behaviour as in Eq. (19), but with 

A^Af = 1.5471 , 2n 2 + n 3 = -0.7424 . 

Another hybrid approximant can be written as 

$\ x ) = [(1 + A x x) n ^ (1 + A 2 x) n * + A 3 x 3 ]" 3 . (20) 

Defining the parameters by means of the accuracy-through-order procedure, we require that 
approximant (20) be real and finite. Typical parameters are 

A l = 0.8724 , A 2 = 0.3360 , A 3 = 0.0197 , 
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rii = 1.7396 , n 2 = 0.3551 , n 3 = -0.3054 
For large x, this yields the same form as in Eq. (16), but with 



A3 3 = 3.3180 , 3n 3 = -0.9162 . 

Approximants f£(x) and fQ j \x) are of comparable accuracy. The existence of multiple 
solutions for the control parameters is somewhat compensated by the mutual closeness of 
the approximants corresponding to different parametric solutions. The accuracy can be 
essentially improved by invoking the minimal difference condition [9]. 



4 Quartic Anharmonic Oscillator 

Let us consider the one-dimensional quartic anharmonic oscillator with the Hamiltonian 

in which x G (—00, +00) and g > 0. The related ground-state energy, obtained by means of 
perturbation theory, is the asymptotic series 

E{g)~Y. a ng n (22) 

n 

in powers of the coupling g. The coefficients a n can be found in Refs. [19,36]. The strong- 
coupling asymptotic behaviour is 

E(g) ~ 0.667986 g 1/3 (g -> 00) . (23) 

The convergence of the factor approximants E% k (g) was shown in [19]. Here we shall compare 
the factor, root, and hybrid approximants with each other. 

We shall consider expansion (22) up to the 6-th order in g. The corresponding factor 
approximant is 

E* 6 (g) = \{1 + A X 9 r (1 + A 2 g) n - (1 + ^3 9T 3 ■ (24) 

The control parameters are uniquely defined from the accuracy-through-order procedure, 
which yields 

A x = 26.4702 , A 2 = 12.4688 , A 3 = 3.8380 , 
ni = 1.8017 x 10~ 3 , n 2 = 0.0547 , n 3 = 0.2005 . 
The strong-coupling limit of Eq. (24) is 

E* 6 (g) ~ i A? A? A? g *+»+* (g - 00) , (25) 

where 

i A^A^Af = 0.7561 , ni +n 2 + n 3 = 0.2570 . 
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Comparing Eqs. (25) and (23), we see that the amplitude in Eq. (25) is predicted with an 
error of 13%, and the power, with an error —23%. 

Dealing with the root and hybrid approximants, with the control parameters defined by 
the accuracy-through-order procedure, we, as always, confront the problem of nonuniqueness 
of solutions. Of course, we shall again reject those solutions which do not guarantee that the 
related approximant be real-valued and finite for finite g. Moreover, we shall present below 
only those root and hybrid approximants of the given order, whose strong-coupling limit is 
closer to the g 1 ^ 3 law. 

Note, first, that the best 5-th order root approximant 



1 3 

2 + 4 



[1+A, g) ni + A 2 if 



n2 



(26) 



with the parameters 

Ax = 24.1009, A 2 = 125.3648, 
has the strong-coupling limit 

E <o) ~ 1 A 2 ' 2 g 2n2+1 



n x = 0.8859 , n 2 = -0.1639 , 



(9 



00 



(27) 



where 



- A 2 2 = 0.3397 , 2n 2 + 1 = 0.6721 



Hence, the amplitude is given with an error of —49% and the power, with an error of 102%. 
The 6-th order root approximant is 



E ( i\g) = =-{[{l + A 1 gr + A 2 gf 2 + A 3 



(28) 



The best accuracy is provided by the parameters 

A 1 = 16.0451 , A 2 = 52.5504 , 

rai = 0.8682 , n 2 = 5.4769 , 
In the strong-coupling limit, this gives 



An 2 n a 2n 2 n 3 

6 ^ 2 y 



A 3 = 37.0388 
n 3 = 0.0197 . 



(29) 



with 



_ A™ 3 = 0.7660 , 
2 2 



2n 2 n 3 = 0.2154 . 



The error of the amplitude is 15% and that of the power is —35%. 
Another solution for the parameters, corresponding to form (28), is 

A 1 = 26.6927 , A 2 = 234.0099 , A 3 = 695.5007 

ni = 0.9638 , n 2 = 0.8934 , n 3 = 0.0653 . 
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This results in the strong-coupling behaviour as 



(9 



oo 



with 



- A T ^ = 0.7664 , 



3n 3 = 0.1958 . 



The errors of the amplitude and power are 15% and —41%, respectively. 
The best hybrid approximant of the form 



7l 2 



+ g) 



nz 



possesses the parameters 

A 1 = 30.9204 , 

ni = 0.9754 , 
The related strong-coupling limit is 



A 2 = 245.4475 
n 2 = 0.0207 , 



Ef (<?)- 2 A ? A ? 9 2n2+n3 



A 3 = 4.1366 
n 3 = 0.2120 . 

(g -»• oo) , 



where 



- A a2 A n " 



0.7570 , 2n 2 + n 3 = 0.2533. 



The amplitude and power errors are 13% and —24%, respectively. 
Another hybrid approximant of the type 



[1 + A 1 9 r (1 + A 2 g) n2 + A 3 g 



"3 



A 2 = 25.9485 , 
n 2 = 0.0263 , 



A 3 = 97.8519 
n 3 = 0.0853 , 



(30) 



(31) 



(32) 



(33) 



Et\ 9 ) = \ 
with the best parameters 

Ax = 7.7952 , 
ni = 2.1670 , 

has the same strong-coupling limit as (30), but with 

i A3 13 = 0.7394 , 3n 3 = 0.2560 . 

The extrapolation of the amplitude is done with an error of 11% and that of the power, with 
an error of —23%. 

In this way, the extrapolation accuracy of the best hybrid approximant (33) is practically 
the same as that of the factor approximant (24). But the latter has the advantage of 
being uniquely defined. The accuracy of factor approximants can be further improved by 
considering either the simple Cesaro averages for the neighbouring approximants, as 



E^ig) + E* k (g) 



(34) 



9 



or the weighted averages of type (8). As we have checked, the accuracy of the weighted 
averages 

E* k (g) = Pk-i(g)E* k _ 1 (g)+Pk(9)E*M > (35) 

with probabilities (9), is essentially better than that of the simple averages (34). One more 
possibility for improving the accuracy is to use the minimal-difference condition for the 
subsequent approximants [9,37]. 



5 Boxed Quantum Particle 

The problem of defining the ground state energy of a quantum particle in a one-dimensional 
box can be formulated [38,39] as the problem of finding the limit of the function 



as g — > oo. The energy is written here in dimensionless units. One has 

E(oo) = ^ = 0.077106 . (37) 
We shall analyse how this value can be extrapolated from the weak-coupling expansion 

E{g)^^- 2 "£a n g n (g ^ 0) . (38) 
The initial coefficients of the latter expansion are 

o = 1 , ai = ^ = 2-467401 , a 2 = ^ = 3.044034 , a 3 = ^ = 1.877713 , 

a 4 = , a 5 = -0.714478 , a 6 = , a 7 = 0.543724 , 
a 8 = , a 9 = -0.517223 , a 10 = , a n = 0.551056 , a 12 = . 
The even factor approximants are 

k 



E ™(9) = ^- 2 na+Asr- (39) 



To guarantee a finite limit, as g — > oo, we impose the restriction 

k 

£n l = 2. (40) 
i=i 

The control parameters Ai and are uniquely defined from the accuracy-through-order 
procedure. 

In the fourth order, we have El(g) with the parameters 

A 1 = 0.30843 + 0.81602 % , A 2 = A\ , m = 1 - 1.13389 % , n 2 =n\. 
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Then we get 

El(oo) = 0.14968 , (41) 

with an error of 90%. 

For the factor approximant E^(g), we find 

A l = 0.44119 , A 2 = 0.08783 + 1.02776 i , A 3 = A\ , 

ni = 1.43469 , n 2 = 0.28266 - 0.86829 i, n 3 = n* 2 . 

As a result, 

E;(oo) = 0.05257 , (42) 

whose error is —32%. 

The approximant E£(g) possesses the parameters 

A 1 = 0.27455 + 0.39441 i , A 2 = A{ , A 3 = 0.03388 + 1.11925 i , A A = A* 3 , 

ni = 0.88578 - 0.77500 i , n 2 = n* , n 3 = 0.11422 - 0.60842 i , n 4 = n* 3 . 
This yields 

^(oo) = 0.10285, (43) 

with an error of 33%. 

For the approximant E* (g), the parameters are 

A 1 = 0.14438 + 0.65311 % , A 2 = A\ , 

A 3 = 0.01626 + 1.16211 i , A 4 = A* 3 , A 5 = 0.29557 , 

ni = 0.30513 - 0.70804 i , n 2 = n* , 
n 3 = 0.06119 - 0.46375 i , n 4 = n* , ra 5 = 1.26735 . 
The limiting value is 

£* (oo) = 0.06201 , (44) 

which deviates from the exact value (37) by an error of —20%. Thus, the sequence of the 
approximants E 2k (oo) displays numerical convergence to the exact result (37). 

Constructing the root approximants, with the control parameters defined by the accuracy- 
through-order procedure, we select, as early, the best approximants. And again, we impose 
the condition that the energy at g — > oo is finite. For example, the root approximant 

E?\g) = [(1 + A 1 9 r + M gf 2 (45) 

is assumed to satisfy the conditions 

n 1 <2, n 2 = l (46) 
guaranteeing that £4^(00) is finite. The parameters of Eq. (45) are 

A 1 = 3.4826 , A 2 = 4.2965 , n x = 0.7085 , n 2 = 1 . 
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The sought limit is 

^ 1} (oo) = 0.0544 , (47) 

which has an error of —30%. 

As an odd root approximant, we consider 



E s\g) = v^r-, \ l +'-r9 (1 + ^1 g) 111 + A 2 g 2 ) , (48) 



1 L 7T 



2 



8ttV 
under the restriction 



ni < 2 , n 2 = - . (49) 



The best solution for the parameters results in 

^(oo) = 0.0648 , (50) 

which has an error of —16%. 

The accuracy of the root approximants is about the same as that of the factor approxi- 
mants. Though the best root approximants are slightly better than the factor approximants 
of the same order, this advantage is spoiled by the problem of multiplicity of solutions for 
the root parameters, while the parameters for the factor approximants are uniquely defined. 

Note that the method of Pade approximants also incorporates the problem of nonuniquely 
defined approximants, since for a given order k of series (2), one may construct a table of 
k Pade-approximants f\M/N], with all M and N satisfying the equality M + N = k. One 
often inclines to use the diagonal Pade approximants [2] , which are not necessarily the most 
accurate. Thus, for series (38) the best Pade approximant of 6-th order is E[ 4 / 2 ](g), yielding 

£ [4/2] (oo) = 0.03855 , 

whose error is —50%. 

The accuracy of the factor approximants can be drastically improved by constructing 
their averages. For instance, even for the simple Cesaro averages we have 

X - (El + E%) = 0.10113 (31%) , 

i (E; + £*) = 0.07771 (0.8%) , 

\ (E* 8 + E* 10 )= 0.08243 (7%), 

where El = E%(oo) and the percentage errors are shown in brackets. Even better is the 
accuracy of the weighted averages (8) involving two nearest neighbours. We find 





= 0.0778 


(0.9%) , 


K 


= 0.0774 


(0.4%) , 


E 10 


= 0.0776 


(0.6%) 



In this way, it looks that the most convenient technique would be by constructing the factor 
approximants and forming their weighted averages. 
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6 Fluctuating Fluid Membrane 

A problem, mathematically very similar to that considered in the previous section, is the 
determination of the pressure of a tensionless membrane between walls [40]. This pressure, 
in dimensionless units, can be presented as the strong-coupling limit 

P(oo)= limP^) (51) 
of a function P(g) that can be calculated for small g — > as a series 



P(g)^^ 2 Y,^g n (</->(>). (52) 



Monte Carlo estimates give [41,42] the value 

P(oo) = 0.0798 ±0.0003 . (53) 

The coefficients in series (52) can be written [39] as 

a = 0.0506606 , a : = 0.125000 , a 2 = 0.154213 , 

a 3 = 0.105998 , a 4 = 0.026569 , a 5 = -0.034229 , a 6 = -0.083251 . 
The even factor approximants are defined by the formula 

P2 k (g) = 2~2 na+^r, (54) 
^9 i=i 

in which, in order to guarantee the finite limit (51), one should set 

k 

E^ = 2 - ( 55 ) 
i=i 

Then one has 

W = 7K- ( 56 ) 
4 i=i 

However, the factor approximants (54), under condition (55), do not provide good accu- 
racy for this problem. For example, P 2 *(oo) = 0.019, which is too small. To the contrary, 
P 4 *(oo) = 0.312 is too large. 

The root approximants here work much better. Thus the approximant 



p 4%) = ^[(i+^<?r+wr > ( 57 ) 

with the parameters 

A 1 = 3.5607 , A 2 = 4.3928 , m = 0.6930 , n 2 = 1 , 
gives P^\oo) = 0.0556, where an error is —30%. 



«2 
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The odd root approximant 



P 5%) = ^{l + Ag [(1 + A 9 r + A 2 gf 2 } , 

where 

A = 2.4674 , Ai = 3.7056 , A 2 = 4.7456 , 

m = 0.6659 , n 2 = 0.5 , 

yields P^\oo) = 0.0681, with an error of —15%. 
The 6-th order root approximant 



(58) 



with 



A ± = 4.8198 , A 2 = 11.4910, A 3 = 13.2536 , 
ni = 0.9893 , n 2 = 0.7762 , n 3 = 2/3 , 

results in P 6 (1) (oo) = 0.0709, which has an error of —11%. 
The odd root approximant 

P$ 1) (9) = ^- 2 (1 + A 9 {{(l + A l9 r + A 29 f>+A 3g z} n3 ) , 



(59) 



(60) 



in which 



A = 2.4674 , Ai = 4.8298 , A 2 = 11.9969 , A 3 = 15.0003 , 

rai = 0.9994 , n 2 = 0.7668 , n 3 = 1/3 , 

gives Py^oo) = 0.07707, which is accurate to within —3.4%. 

Recall that the root approximants are not uniquely defined from the accuracy-through- 
order procedure. Here only the best of them are presented. Different solutions for the control 
parameters usually lead to the approximants that are close to each other. For instance, 
another solution for form (60) would give P 7 (g) with the parameters 

A = 2.4674 , Ai = 2.3970 , A 2 = 6.1342 , A 3 = 14.9918 , 

ni = 2.0166 , n 2 = 0.7657 , n 3 = 1/3 . 

This gives practically the same limit P 7 (2) (oo) = 0.07706 as P 7 (1) (oo). 

Pade approximants, invoked for this problem, either contain divergencies, or yield prin- 
cipally incorrect results with the negative values of pressure. 
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7 One-Dimensional Antiferromagnet 

The ground-state energy of the one- dimensional Heisenberg antiferromagnet is known exactly 
[43], 

E = -0.4431 . (61) 
From another side, this energy can be considered as the limit 

E = lim E(t) 

of the temporal energy E(t). The latter can be calculated for small t, yielding the so-called 
t-expansion [44] 

E(t)c±-\j2*nt n (f->0), (62) 

in which 

16 

a = 1 , ai = 4 , a 2 = -8 , a 3 = - y , a 4 = 64 . 
The even 4-order factor approximant extrapolating expansion (62) is 

Et{t) = -l(l + A 1 (1 + A 2 tr , (63) 

which in the long-time limit gives 

E* 4 (t) ~ - i A\ X A 2 2 t m+m (i -> oo) . 

The parameters of Eq. (63) are uniquely defined by the accuracy-through-order procedure, 
with the condition 

^i + n 2 = , 

guaranteeing the finiteness of E2(oo). This gives ^4(00) = —0.570, with an error of 29%, as 
compared to the Hulthen result [43]. 
The odd factor approximant 

E* 5 (t) = - i [1 + 4i(l + Ax t) ni (l + A 2 t) n >] (64) 

has the limiting behaviour 

E* 5 (t) ~ - J (l + 4A™M™ 2 r i+ " 2+1 ) (t -> 00) , 

with the restriction 

™i + n 2 + 1 = . 

From here, E£(oo) = —0.4452, which is accurate to within 0.5%. 
The best root approximant 

Eg\t) = -\ {l + 4t[(l + A 1 t) ni + A 2 1?] r>a } , (65) 
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with the asymptotic form 

Ei j) (t) ~ - ^ (l + 4A™ 2 t 2n2+1 ) (i -> oo) , 

requires the condition 

2n 2 + 1 = . 

The related limit is E§\oo) = —0.4743, whose error is 7%. 

The diagonal Pade approximant E[ 2 /2](t) gives E[ 2 /2](oo) = —0.3289, with an error of 
26%. 

8 Combined Factor-Exponential Approximants 

When the asymptotic behaviour of a function at large x — > oo is a combination of the power- 
law and exponential dependence, the extrapolation of the corresponding series at small x — > 
can be done by combining the factor and exponential approximants. We shall illustrate this 
by considering the Airy function satisfying the Airy equation 

d 2 

— Ai(x) - xAi(x) = . (66) 
ax 2 

The solution to this equation at small x can be written as a series 

Ai{x) ~£a„ * n ( x ^ °) > ( 67 ) 

n 

substituting which in Eq. (66) gives the coefficients 

a = 0.355028 , a x = -0.258819 , a 2 = , a 3 = — , 

6 

ai a ai 

12 ' 180 ' 504 ' 8 

At large x, the Airy function behaves as 

Ai(x) ~ x- 1/4 exp f- ^ x 3/2 ") (x -> oo) . (68) 
2^/71 V 3 / 

The combined factor-exponential approximant of fifth order can be represented as 

A* 5 (x) = ao(l + a;)" 1 exp {bx 2 (l + A 2 x) ri2 } , (69) 

where the methods of Refs. [11,12] are involved. The parameters, found from the accuracy- 
through-order procedure, are 

Ai = 1.480028, A 2 = 1.400808, b = -0.805208 , m = -0.492565 , n 2 = -0.505184 . 

The large-x behaviour of approximant (69) is 

A* 5 (x) ~ aoA™ 1 x ni exp (bA^ x n2+2 ) , (70) 

where 

bA r 2 12 = -0.697141 , n 2 + 2 = 1.494816 , 

which are very close to the values —2/3 and 3/2, respectively. Approximant (69) represents 
the exact Airy function very well. Thus, the error for x < 10 is less than 1%. 
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9 Discussion of Other Problems 



The considered technique can be applied to any problem requiring an extrapolation of small- 
variable asymptotic series to the large-variable region. We have accomplished such an ex- 
trapolation for a variety of problems. However, we would not like to overload this paper by 
the description of calculational details corresponding to the related problems. Instead, we 
shall just briefly summarize the cases we have analysed. 

The luminescent intensity of donor-acceptor recombination [45] was extrapolated by 
means of the factor approximants in Ref. [19]. Now, we have also considered the extrapola- 
tion with the help of the root approximants. The best variants of the latter are not better 
than the factor approximants of the same order. The accuracy can be improved by invok- 
ing averages (8). The best Pade approximants are much less accurate than the self-similar 
approximants. 

The nonlinear Schrodinger equation is met in several physical applications. A very impor- 
tant and interesting application is the description of coherent fields of trapped Bose atoms, 
when the corresponding nonlinear equation is called the Gross-Pitaevskii equation [46-50]. 
Self-similar approximants for the ground-state wave function and energy of the latter equa- 
tion were considered in Refs. [16,51] and for the whole spectrum of excited energies in Refs. 
[52,53]. Here, we have compared the extrapolation procedure, based on both factor and 
root approximants, for the ground-state energy. We find that for this problem the factor 
approximants are more accurate than the root approximants. 

The factor approximants are also convenient for describing critical phenomena. Ther- 
modynamic characteristics in the vicinity of the critical point exhibit the behaviour typical 
of one of the factors [54,55]. Earlier, we checked the applicability of subsequent factor 
approximants to different critical phenomena [18,19]. Now we have shown that the hy- 
brid factor-root approximants can also be used for describing critical phenomena. We have 
checked this for the elliptic integral with logarithmic singularity [19] and for the so-called 
(2+1) dimensional Ising model [56]. The accuracy of the approximants can be improved by 
involving the minimal- difference condition [9,37]. 

Critical behaviour may also appear in the solutions to nonlinear differential equations 
[23,57-59]. We studied this effect for the Ruina-Dietrich equation by using the factor ap- 
proximants [18]. As we have now checked this critical behaviour can also be described by 
the hybrid factor-root approximants. 

In the majority of the cases we have investigated, Pade approximants are essentially less 
accurate than the self-similar approximants. Often, Pade approximants are not applicable 
at all, being divergent or qualitatively incorrect. In particular cases, it is possible to fit the 
correct behaviour of a sought function by making manipulations with Pade approximants 
raising them to some fractional powers [60] . However, such a way of fitting is too arbitrary to 
serve as a serious method. In addition, this fitting actually results in constructions analogous 
to self-similar root approximants, though slightly spoiled and not so symmetric. 
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10 Conclusion 



A power series in powers of an asymptotically small variable x — > can be effectively 
extrapolated to the region, where this variable is large, and even to the limiting case x — > 
oo. This can be done by means of the self-similar factor and root approximants. Hybrid 
approximants, combining factors and roots, can also be used for extrapolation. The control 
parameters can be defined from the accuracy-through-order relations. This procedure yields 
unique solutions for the factor approximants but multiple solutions for root approximants. 
Fortunately, the multiplicity of solutions for the root parameters is often not as dangerous, 
leading to approximants that are very close to each other. 

The extrapolation by the factor approximants is preferable, being more accurate, when 
the sought function either increases to infinity, as x — > oo, or diminishes to zero. When the 
function tends, as x — > oo, to a nonzero finite value, the accuracy of the root approximants 
can become higher than that of the factor approximants. 

The usage of the complete root-factor technique could be justified in those cases when 
the number of available terms in a small-variable asymptotic expansion is limited, if the 
derivation of the higher-order terms is too costly or even impossible at all. The trouble with 
multiple solutions could be overcome by imposing additional restrictions on the behaviour 
of solutions at infinity. 

A very important feature of the root and hybrid root-factor approximants is their non- 
trivial behaviour at infinity allowing for defining the so-called corrections to scaling, while 
the direct usage of the factor approximants yields only the leading exponent. This problem 
will be studied in detail in a separate paper. 
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